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Abstract. We summarize recent theoretical developments in the field of radiative and semileptonic penguin decays. 
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AMBIGUITY OF THE NEW PHYSICS SCALE 

Within the indirect search for New Physics (NP) there is a an ambiguity of the new physics scale. In the model- 
independent approach using the effective electroweak hamiltonean, the contribution to one specific operator ff\ can be 
parametrized via 

(C; M /M w + C' NP /A NP ) x e% (1) 

where the first term represents the SM contribution at the electroweak scale Mw and the second one the NP contribution 
with an unknown coupling C' SM and an unknown NP scale Anp- 

The radiative and semi-leptonic penguin modes, b — > sy and b — > s£ + £~ , are flavour changing neutral current (FCNC) 
processes and, thus, are highly sensitive for new degrees of freedom via virtual effects (for reviews, see Ull^yl]) The 
non-existence of large NP effects in flavour observables in general fl.Ht] implies the famous flavour problem, namely 
why FCNC are suppressed. Either the mass scale of the new degrees of freedom, Anp, is very high or the new flavour- 
violating couplings, C' NP , are small for (symmetry?) reasons that remain to be found. For example, assuming generic 
new flavour-violating couplings of 0( 1 ), the present data on K-K mixing implies a very high NP scale of order 10 3 -10 4 
TeV depending on whether the new contributions enter at loop- or at tree-level. In contrast, theoretical considerations 
on the Higgs sector, which is responsible for the mass generation of the fundamental particles in the SM, call for NP at 
order 1 TeV. As a consequence, any NP below the 1 TeV scale must have a non-generic flavour structure. In addition, 
also the present electroweak data indicate a slightly higher NP scale, a data-driven problem known as little hierarchy 
problem. 

The present measurements of B decays, especially of FCNC processes, already significantly restrict the parameter 
space of NP models. In general such bounds from flavour physics are model-dependent, but often much stronger than 
the ones derived from other measurements. In any case, the indirect flavour information will be most valuable when 
the general nature and the mass scale of NP will be known. 

THE INCLUSIVE DECAY B -> X S)d y 

Perturbative contributions: Among the rare decay modes, the inclusive decay B — > X s y is the most important one, 
because it is theoretically well-understood and at the same time it has been measured extensively at the B factories. 
While non-perturbative corrections to this decay mode are subleading and recently estimated to be well below 
10% [6], perturbative QCD corrections are the most important corrections. Within a global effort, a perturbative QCD 
calculation to the next-to-next-to-leading-logarithmic order level (NNLL) has quite recently been performed and has 
led to the first NNLL prediction of the B — > X s y branching fraction [7] with a photon cut at Ey = 1.6GeV (including 
the error due to nonperturbative corrections): 

£§(B^X s y) NNLL = (3. 15 ±0.23) x 1(T 4 . (2) 
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This result is based on various highly-nontrivialperturbative calculations jl||,|lMlHEQBB[I3,[Ii- The 
combined experimental data leads to (Heavy Flavor Averaging Group (HFAG) [19]) 

■S§{B^X s y) = (3.55 ±0.24 ±0.09) x 10~ 4 , (3) 

where the first error is combined statistical and systematic, and the second is due to the extrapolation in the photon 
energy. Thus, the SM prediction and the experimental average are consistent at the 1.2(7 level. This is one important 
example that the CKM theory is not only confirmed by the data entering into the CKM unitarity fit, but also by many 
additional flavour mixing phenomena. 

Nonperturbative contributions: It was noted long ago [20], that there is no local OPE for the inclusive decay B — > X s y 
if one considers operators beyond the leading electromagnetic dipole operator Then, there are so-called resolved 
photon contributions which contain subprocesses in which the photon couples to light partons instead of connecting 
directly to the effective weak-interaction vertex [6]. Only recently, a systematic analysis 16J] of all resolved photon 
contributions related to other operators in the weak Hamiltonian has established this breakdown of the local OPE 
within the hadronic power corrections as a generic result. Clearly, estimating such nonlocal matrix elements is very 
difficult, and an irreducible theoretical uncertainty of ±(4 — 5)% for the total CP averaged decay rate, defined with 
a photon-energy cut of Ey = 1.6 GeV, remains [6]. This result strongly indicates that the theoretical efforts for the 
B — > X s y mode have reached the nonperturbative boundaries, but it also reconfirms the dominance of the perturbative 
contributions. 

There is another positive result induced by this new analysis: Until recently, it was believed that the long-distance 
contributions from the intermediate w-quark in the penguin loops are critical. They are suppressed in the B — > X s y mode 
by the CKM matrix elements. InB^ X^y, there is no CKM suppression, and one must account for the nonperturbative 
contributions that arise from the M-quark loops. A simple dimensional estimate leads to an uncertainty of at least 10%. 
However, this contribution vanishes in the total CP-averaged rate of B — > X s y at order A /my [6]. This result applies 
to the total rate of B — >• Xjj as well. Thus, there is no power correction due to the M-quark loops in the total rate of 
B — > X £ /Y at order A/nth, which implies that the CP-averaged decay rate of B — > X c [j is as theoretically clean as the 
decay rate of B — > X s y U. The present NLL-prediction reads lf2~ll l22tl (for Ey > 1 .6 GeV): 

3S{B^X d y) = (\.38 i^|^±0.15 C KM±0.09 param .±0.05 SC aieJ X lO" 5 . (4) 

The large CKM-uncertainty is due to V tc i. The uncertainty due to the charm-scheme dependence can be reduced 
significantly by a NNLL analysis in analogy to B — > X s y. Already the NLL prediction in combination with the first real 
measurement of this inclusive mode l23ll leads to interesting bounds on NP lf24ll . 

But there is also a negative consequence: Without CP-averaging the non-perturbative the M-quark contribution 
survives and dominates the direct CP asymmetry inB^ X s( /y. Perturbative contributions exhibit a triple suppression; 
one has an a s suppression in order to have a strong phase, a Cabibbo suppression of A 2 , and a GIM suppression 
of (m c /my) 2 A/my reflecting the fact that in the limit m c —> the direct CP violation vanishes. Thus, the non- 
perturbative (resolved) M-quark contribution is enhanced by a factor 1 / a s compared to the perturbative contributions. 
Model estimates of the resolved contribution lead to l25ll 

- 0.6% < stf(B ->■ X s y) SM < 2.8%, (5) 

which covers most of the experimentally allowed range fl^l : &/(B — > X s y) = —(1.2 ±2.8)% . But the untagged direct 
CP asymmetry survives. All resolved contributions cancel at order A/my, so we still have the zero-prediction for the 
untagged CP asymmetry for B — > X s+d y l26l l27i l2lll . This prediction is directly based on the CKM-unitarity and on 
the smallness of the U-spin breaking parameters (m s /m;,) 2 . Thus, it represents a clean test whether new CP phases are 
active or not. 

There are still some open issues in the decay B — > X s y which ask for further study. First, there are three-loop 
matrix elements of the leading four-quark operators, which have first been calculated within the so-called large-/3o 
approximation lfl7l I28H29I1 . A calculation that goes beyond this approximation by empl oying a n interpolation in the 
charm quark mass m c from m c > my to the physical m c value has been presented in Refs. IU8l 13011 . In this interpolation 
the a 2 /3o result lfl7ll is assumed to be a good approximation for the complete a 2 result for vanishing charm mass. 
It is this part of the NNLL calculation which is still open for improvement. Indeed, there are several collaborations 
presently working on this issue. A complete calculation of the three-loop matrix elements of the four-quark operators 
ffi 2 for vanishing charm mass is in progress and will cross-check the error estimate due to the interpolation Oil l32ll . 



Partial results are already available l33ll . In addition, there is an effort l34ll to calculate the matrix elements for arbitrary 
m c directly. All these efforts will finally eliminate the 3% uncertainty due to the interpolation within the present NNLL 
prediction Q]. 

Second, in the measurement of the inclusive mode B —> X> y one needs cuts in the photon energy spectrum to suppress 
the background from other B decays which induces additional sensitivities to non-perturbative physics. These shape- 
function effects were taken into account in the experimental analysis. The corresponding theoretical uncertainties 
due to this model dependence are reflected in the extrapolation error of the experimental results. The extrapolation 
is done from the experimental energy cut values down to 1.6 GeV. But a cut around 1.6GeV might not guarantee 
that a theoretical description in terms of a local OPE is sufficient because of the sensitivity to the additional scale 
A = mb — 2Ey l35ll . A multiscale OPE with three short-distance scales wj, y/m^A, and A has been proposed to connect 
the shape function and the local OPE region. Recently, such additional perturbative cutoff-related effects have been 
calculated to NNLL precision by the use of SCET methods |36ll37ll38l1 . Such perturbative effects due to the additional 
scale are negligible at l.OGeV but of order 3% at 1.6 GeV 13611 . The size of these effects at 1.6 GeV is similar to the 3% 
higher-order uncertainty in the present NNLL prediction. However, the numerical consistency of the SCET analysis has 
recently been questioned 13911 . Far away from the endpoint (£0 = 1 .6 GeV), the logarithmic and nonlogarithmic terms 
cancel; this kind of cancellation was already observed in Ref. [40]. Within the resummation of the cutoff-enhanced 
logarithms this feature leads to an overestimate of the 0(0$) terms 1I39I1 . Further work is needed to clarify this issue. 

Third, the B —> X s y decay rate is normalized to the charmless semileptonic rate in order to separate the charm 
dependence. Then the quantity C — \V u b\ 2 / \V c b\ 2 x T[B — > X c ev]/T[B — > X u ev] enters the theoretical prediction. 
Recently, a scheme dependence in the determination of the prefactor C was noticed 141], which is around 3%, thus 
within the perturbative uncertainty |39]. The two determinations in the IS scheme [42] and in the kinetic scheme i4lll 
differ through renormalization schemes, methodology, and experimental input. An update of the analysis within the 
IS scheme might resolve part of the corresponding uncertainty. 

THE INCLUSIVE DECAY B -> X s £ + £~ 

Perturbative contributions: This inclusive mode is also dominated by perturbative contributions, if one eliminates cc 
resonances with the help of kinematic cuts. In the so-called 'perturbative g 2 -windows' below and above the resonances, 
namely in the low-dilepton-mass region 1 GeV 2 < q 2 = m 2 e < 6 GeV 2 , and also in the high-dilepton-mass region 
with q 2 > 14.4 GeV 2 , theoretical predictions for the invariant mass spectrum are dominated by t he p erturbative 
contributions. QCD corrections are calculated to NNLL precision JH |13,|4l|4^|43|4il|4§|5^|5ll|52]. More 
recently electromagnetic corrections were calculated: NLL quantum electrodynamics (QED) two-loop corrections 
to the Wilson coefficients are of 0(2%) l49ll . Also, in the QED one-loop corrections to matrix elements, large 
collinear logarithms of the form \og{m\/m 2 ) survive integration if only a restricted part of the dilepton mass spectrum 
is considered. These collinear logarithms add another contribution of order +2% in the low-g 2 region for B —> 
X SJ u + jU~ JH]. For the high-g 2 region, one finds -8% JH]. 

Quark-Hadron-Duality : The integrated branching fraction of B —> X s £ + l~~ is dominated by the resonance back- 
ground which exceeds the nonresonant charm-loop contribution by two orders of magnitude. As has been recently 
noticed 115411 . this feature should not be misinterpreted as a striking failure of global parton-hadron duality which pos- 
tulates that the sum over the hadronic final states, including resonances, should be well approximated by a quark-level 
calculation [55]. Crucially, the charm-resonance contributions to the decay B — > X s £ + £~ are expressed in terms of 
a phase-space integral over the absolute square of a correlator. For such a quantity global quark-hadron duality is 
not expected to hold. Nevertheless, local quark-hadron duality (which, of course, also implies global duality) may be 
reestablished by resumming Coulomb-like interactions 15411 . 

THE EXCLUSIVE DECAY B -»■ K^i+T 

QCD factorization: The method of QCD factorization (QCDF) JH] and its field-theoretical formulation SCET ll57Ll58ll 
form the basis of the up-to-date predictions of exclusive B decays. There is also a factorization formula for the exclusive 
semileptonic B decays, such as B — > K*£ + £~ [59, 60]. The hadronic form factors can be expanded in the small ratios 
A/ m/, and A/E, where E is the energy of the light meson. If we neglect corrections of order l/nij, and a s , the seven 
a priori independent B — > K* form factors reduce to two universal form factors t;± and &i 116 ill . These theoretical 



simplifications of the QCDF/SCET approach are valid to the kinematic region in which the energy of the K* is of the 
order of the heavy quark mass, that is, q 2 <C m 2 B . Thus, factorization formula applies well in the dilepton mass range 
1 GeV 2 < q 2 < 6 GeV 2 . 

Full angular analysis: In the near future, a full angular analysis will become possible. This rich information 
allows for the design of observables with specific NP sensitivity and reduced hadronic uncertainties ll62il63ll . These 
observables are constructed in such a way that the soft form factor dependence cancels out at leading order in a s and in 
A / ntb for all low dilepton masses, and they have much higher sensitivity to new right-handed currents than observables 
that are already accessible via the projection fits 1164 l62l 16311 . In these optimized observables, the unknown A/m& 
corrections are the source of the largest uncertainty. Further detailed NP analyses of such angular observables have 
been presented in Refs. Ir55l l66ll . A full angular analysis provides high sensitivity to various Wilson coefficients, but 
the sensitivity to new weak phases is restricted, mainly due to large experimental uncertainties 1671 16311 . Observables 
defined at high-g 2 also represent very interesting observables due to the fact that A/nif, corrections can be estimated 
with the help of heavy-quark effective theory (HQET), but formfactors at high-g 2 have to be presently extrapolated 
from the \ow-q 2 region jr5^ |69l |70ll . 
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